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"Aerodynamic  Flow  Field  Approximations,"  and  24040775,  "Lifting  Entry  Config¬ 


urations."  The  study  period  was  February  1983  to  June  1986. 


The  experimental  program  described  in  this  report  produced  a  very  large 
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SECTION  1.0 


INTRODUCTION 

The  prediction  of  the  aerodynamic  characteristics  of  the  latest  missile 
configurations  being  studied  has  involved  elliptical  cross-section  missile 
bodies.  Past  efforts  in  this  area  have  shown  deficiencies  in  predicting  the 
aerodynamic  characteristics  of  these  type  of  configurations.  The  first  task 
of  an  AFWAL/FIMG  contracted  effort  entitled,  "Aerodynamic  Analysis  for  Mis¬ 
siles"  was  the  evaluation  of  10  aerodynamic  prediction  methods  for  four 
classes  of  missile  configurations.  The  lifting  missile  class  consisted  of 
elliptical  body  configurations  with  wings  and  tails.  The  limited  comparisons 
made  of  the  Supersonic-Hypersonic  Arbitrary  Body  Program  (S/HABP)  with  the 
elliptical  bodies  for  Mach  numbers  from  2.0  to  4.0  showed  poor  results, 
particularly  at  the  lower  Mach  numbers. 

An  in-house  work  element  was  then  initiated  to  more  completely  determine 
which  methods  or  combination  of  methods  available  in  the  S/HABP  code  could 
give  acceptable  results  for  this  type  of  missile  body.  Test  data  for  several 
missile  bodies  ranging  from  circular  to  a  3-to-l  ellipticity  ratio  for  Mach 
numbers  1.5  to  4.63  were  compared  with  the  results  from  the  S/HABP  code  usiny 
a  variety  of  pressure  methods.  The  results  of  the  effort  showed  that  no 
typical  application  of  any  of  the  pressure  methods  in  S/HABP  would  provide 
good  results  across  the  Mach  number  range  and  that  parametric  wind  tunnel 
data,  particularly  pressure  data,  would  be  required  to  determine  the  cause  of 
the  mismatch  of  theory  versus  test. 

To  provide  these  data,  a  series  of  wind  tunnel  tests  were  conducted  in 
the  AEDC  VKF  Tunnel  A  facility  on  basic  elliptical  missile  bodies.  Three 
elliptical  body  models  with  ellipticity  ratios  of  3.0:1,  2.5:1,  and  2.0:1  were 
built  and  tested  at  Mach  numbers  from  1.5  to  5.0.  Both  force  and  moment  and 
pressure  data  were  obtained  as  well  as  flow  visual ization  data  such  as  vapor 
screens,  oil  flow,  and  shadowgraphs. 
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To  complete  the  test  data  base,  the  three  elliptical  body  models  were 
tested  in  the  Aerodynamic  Wind  Tunnel  (4T)  of  the  Propulsion  Wind  Tunnel 
Facility  at  Mach  numbers  from  0.4  to  1.3.  Both  force  and  moment  and  pressure 
data  were  obtained.  In  addition  to  model  surface  pressures,  data  were  also 
obtained  on  two  static  pressure  pipes  mounted  near  the  tunnel  top  and  bottom 
walls  for  one  of  the  configurations  at  Mach  numbers  up  to  1.05.  The  pipe 
pressure  data  were  used  to  determine  if  there  were  any  significant  wall 
interference  effects  on  the  model  surface  pressure  distribution. 

This  report  summarizes  the  results  presented  previously  in  five  AFWAL 
Technical  Memorandums  (References  1  through  5).  The  information  in  the 
sections  on  apparatus  and  test  description  has  been  extracted  from  the  AEOC 
test  reports  (References  6  and  7). 


SECTION  2.0 


APPARATUS 


2.0  Test  Facilities 

AEDC  VKF  Tunnel  A  (Figure  1)  is  a  continuous,  closed-circuit,  variable 
density  wind  tunnel  with  an  automatically  driven  flexible-plate-type  nozzle 
and  a  40-  by  40-inch  test  section.  The  tunnel  can  be  operated  at  Mach  numbers 
from  1.5  to  5.5  at  maximum  stagnation  pressure  from  29  to  195  psia,  respec¬ 
tively,  and  stagnation  temperatures  up  to  750°R  at  Mach  number  5.5.  Minimum 
operating  pressures  range  from  about  one-tenth  to  one-twentieth  of  the  maximum 
at  each  Mach  number.  The  tunnel  is  equipped  with  a  model  injection  system 
which  allows  removal  of  the  model  from  the  test  section  while  the  tunnel 
remains  in  operation.  A  description  of  the  tunnel  and  airflow  calibration 
information  may  be  found  in  Reference  8. 

The  AEDC  Aerodynamic  Wind  Tunnel  (4T)  is  a  closed-loop  continuous  flow, 
variable-density  tunnel  in  which  the  Mach  number  can  be  varied  from  0.1  to  1.3 
and  can  be  set  at  discrete  Mach  numbers  of  1.6  and  2.0  by  placing  nozzle 
inserts  over  the  permanent  sonic  nozzle.  At  all  Mach  numbers,  the  stagnation 
pressure  can  be  varied  from  300  to  3,400  psfa.  The  test  section  is  4-feet 
square  and  12.5  feet  long  with  perforated,  variable  porosity  (0.5-  to  10- 
percent  open)  walls.  It  is  completely  enclosed  in  a  plenum  chamber  from  which 
air  can  be  evacuated,  allowing  part  of  the  tunnel  airflow  to  be  removed 
through  the  perforated  walls  of  the  test  section.  The  model  support  system 
consists  of  a  sector  and  sting  attachment  which  has  a  pitch  angle  capability 
of  -8  to  27  degrees  with  respect  to  the  tunnel  centerline  and  a  roll  capabil¬ 
ity  of  f100  degrees  about  the  sting  centerline.  A  more  complete  description 
of  the  tunnel  may  be  found  in  Reference  8. 

2.2  Test  Articles 

The  test  articles  were  elliptic  missile  body  configurations  with  elli- 
pticity  ratios  of  2.0,  2.5,  and  3.0  to  1.0.  The  three  models  were  designed 
and  fabricated  from  aluminum  at  AEDC,  based  upon  criteria  provided  by 
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AFWAL/FIMG.  The  models  were  power-law  bodies  with  an  exponent  of  0.5  and  had 
the  same  longitudinal  distribution  of  cross-sectional  area.  These  models  were 
based  on  one  of  a  series  of  related  bodies  with  cross-sectional  ellipticity 
tested  by  NASA  at  Langley  Research  Center  (Reference  9).  The  semimajor  and 
semiminor  axis  ordinates  were  derived  from  the  following  equations: 

For  horizontal  projection  (semimajor  axis) 

a  =  amax  *  X0*5 
LM0*5 

and  for  vertical  projection  (semiminor  axis) 
b  =  bmax  *  x0*5 
LM6^5 

Details  of  the  models  are  given  in  Figure  2.  Model  configuration  designation 
is  presented  in  Table  1  and  a  listing  of  the  configurations  tested  is  given  in 
Table  2. 

For  the  pressure  phase,  each  model  was  instrumented  with  191,  0.045  inch 
diameter  surface  pressure  orifices  and  one  base  pressure  orifice.  The  lo¬ 
cation  and  designation  of  the  pressure  orifices  are  presented  in  Figure  3a  and 
were  identical  for  all  three  models.  Table  3  provides  nominal  axial  locations 
from  the  nosetip  and  nominal  radial  locations  from  the  top  ray  (positive 
clockwise  looking  upstream)  of  the  pressure  orifices.  Base  pressure  was 
measured  with  an  orifice  located  halfway  between  the  sting  and  model  outer 
surface  as  shown  in  Figure  3c.  For  the  transonic  pressure  test,  static  pres¬ 
sure  pipes  were  mounted  on  the  centerline  of  the  top  and  bottom  walls. 

Details  of  the  pressure  pipes  are  shown  in  Figure  4.  Each  pipe  had  3D  pres¬ 
sure  orifices  on  each  of  the  model  and  wall  sides  of  the  pipes.  A  more 
detailed  discussion  Is  given  In  Reference  5.  Two  thermocouples  were  attached 
to  the  inner  wall  of  each  model  and  the  general  location  is  given  in  Table  3. 


For  the  Tunnel  A  force  phase,  the  pressure  tubes  were  cut,  sealed  and 
secured  as  required  to  prevent  interference  on  the  balance  measurements.  Base 
pressure  was  measured  with  a  fast-response  pressure  transducer  module  located 
in  a  sting  component  approximately  18  inches  downstream  of  the  model  base.  An 
8-degree  prebend  installation  arrangement  was  used  to  provide  the  angle- 
of-attack  range  from  -4  to  20  degrees.  One  set  of  sting  components  was  used 
for  the  pressure  and  oil-flow  phases  and  a  different,  but  similar  set  of 
components  was  used  for  the  force  and  vapor-screen  phases.  Sketches  of  the 
installation  arrangements  for  both  Tunnel  A  and  4T  are  shown  in  Figure  5. 

2.3  Test  Instrumentation 

During  the  pressure  phases  of  testing,  model  surface  and  base  pressures 
were  measured  with  Pressure  Systems  Incorporated  electronically-scanned 
pressure  modules  referenced  to  a  near  vacuum.  Six  Model  ESP-32  modules  were 
used  for  the  VKF  supersonic  testing;  four  2.5-psid  range  and  two  5-psid  range. 
Each  ESP-32  module  has  32  pressure  ports  with  a  silicon  pressure  transducer 
for  each  port  that  can  be  digitally  addressed  and  calibrated  on-line.  The 
relatively  small  size  of  each  pressure  module  (1.0  X  2.0  X  2.5  inch)  permitted 
on-board  mounting,  which  resulted  in  a  significant  reduction  in  pressure 
stabilization  time  and  a  significant  increase  in  the  data  acquisition  rate. 
Nine  Model  ESP-48  modules,  each  with  a  15-psid  range,  were  used  in  the  4T 
transonic  testing;  five  48-port  modules  located  in  the  model  and  two  48-port 
and  two  16-port  modules  located  outside  the  test  section  and  connected  to  the 
static  pressure  pipes.  Model  wall  temperatures  were  measured  with  two 
Chromel-Alumel  thermocouples  attached  to  the  inner  wall  of  the  model. 

During  the  Tunnel  A  force  phase,  model  base  pressures  were  measured  with 
a  miniature  pressure  transducer  module  manufactured  by  the  Scanivalve  Corpo¬ 
ration.  The  module  contains  eight  fast-response  pressure  transducers  with  a 
range  of  1  psid  referenced  to  a  near  vacuum.  The  transducers  are  diffused 
silicon-diaphragm- type  strain  gage  sensors  fabricated  by  Druck  Incorporated. 
The  small  size  of  the  pressure  module  (0.5  X  1.2  X  1.3  inch)  permitted 
mounting  in  a  sting  component  downstream  of  the  model  base. 


In  Tunnel  A  model  shadowgraph  or  schlieren  photographs  were  obtained  with 
a  double-pass  optical  flow-visualization  system  with  a  35-in-diameter  field  of 
view.  Vapor-screen  photographs  were  obtained  with  two  Hasselbald  70-mm  still 
cameras  mounted  on  the  operating  side  of  the  test  section  and  a  D.B.  Milliken- 
55,  16-mm  movie  camera  at  12  frames  per  second  mounted  on  the  non-operating 
side  of  the  test  section.  Model  flow  field  illumination  was  provided  with  a 
15-mw  helium-neon  continuous  wave  laser  with  the  output  beam  expanded  In  one 
direction  with  a  8-mm  cylindrical  lens.  The  cylindrical  lens  was  rotated  so 
that  the  light  plane  was  perpendicular  to  the  model  centerline  for  each  model 
angle  of  attack. 

Oil -flow  photographs  were  obtained  with  three  Varitron  70-mm  still 
cameras  mounted  on  the  operating  side  of  the  test  section  to  facilitate 
simultaneous  photo  acquisition  along  the  full  length  of  the  model.  An  auto¬ 
matic  camera  control  system  was  used  to  provide  automatic  shutter  sequencing 
at  4-second  intervals.  For  the  second  phase  of  oil  flow  testing,  Hasselbald 
70-mm  still  cameras  were  used  to  photograph  the  models  in  the  access  tank. 
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TEST  DESCRIPTION 


3.1  Test  Conditions 


A  complete  listing  of  test  conditions,  configurations,  and  run  numbers  is 
presented  in  Table  2. 

3.2  Test  Procedures 

3.2.1  General 

In  the  VKF  wind  tunnels,  (A,  B,  C),  the  model  is  mounted  on  a 
sting  support  mechanism  in  an  installation  tank  directly  underneath  the  tunnel 
test  section.  The  tank  is  separated  from  the  tunnel  by  a  pair  of  fairing 
doors  and  a  safety  door.  When  closed,  the  fairing  doors,  except  for  a  slot 
for  the  pitch  sector,  cover  the  opening  to  the  tank  and  the  safety  door  seals 
the  tunnel  from  the  tank  area.  After  the  model  is  prepared  for  a  data  run, 
the  personnel  access  door  to  the  installation  tank  is  closed,  the  tank  is 
vented  to  the  tunnel  flow,  the  safety  and  fairing  doors  are  opened,  the  model 
is  injected  into  the  alrstream,  and  the  fairing  doors  are  closed.  After  the 
data  are  obtained,  the  model  is  retracted  into  the  tank  and  the  sequence  is 
reversed  with  the  tank  being  vented  to  atmosphere  to  allow  access  to  the  model 
in  preparation  for  the  next  run.  The  sequence  is  repeated  for  each  configura¬ 
tion  change. 

The  tunnel  4T  test  section  (4  feet  square  by  12  feet  long)  is 
accessed  through  a  removable  side  wall.  Model  changes  require  shutting  down 
the  tunnel  and  opening  the  side  wall.  For  this  reason,  model  changes  are  done 
after  all  of  the  required  Mach  numbers  are  completed.  For  each  Mach  number 
the  tunnel  conditions  are  held  constant  while  varying  model  attitude.  The 
data  are  recorded  at  selected  angles  using  the  pi tch/rol 1 -pause  technique. 

This  sequence  is  repeated  for  each  configuration. 
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3.2.2  Data  Acquisition  and  Reduction 


3.2.2. 1  Tunnel  A 


Model  attitude  positioning  and  data  recording  were 
accomplished  with  the  point-pause  and  sweep  modes  of  operation,  using  the 
Model  Attitude  Control  System  (MACS).  Model  pitch  and  roll  requirements  were 
entered  Into  the  controlling  computer  prior  to  the  test.  Model  positioning 
and  data  recording  operations  were  performed  automatically  during  the  test  by 
selecting  the  list  of  desired  model  attitudes  and  initiating  the  system. 


Point-pause  force  data  were  obtained  for  finite  values  of  ALPHA  and  BETA 
with  a  delay  before  each  point  to  allow  the  base  pressure  to  stabilize.  Each 
data  point  for  this  mode  of  operation  Is  the  result  of  a  Kaiser-Bessel  digital 
filter  utilizing  16  samples  over  a  time  span  of  0.33  seconds. 


The  continuous  sweep  force  data  were  obtained  for  a  fixed  value  of  PHi 
with  a  sweep  (ALPHA)  rate  of  0.5  deg/sec.  A  data  sample  was  recorded  every 
0.0208  seconds  and  a  Kaiser-Bessel  digital  filter  was  applied  to  every  16 
samples  to  produce  a  sample  data  point  every  0.01  degrees  In  pitch.  The  data 
were  then  Interpolated  to  obtain  the  data  at  the  requested  model  attitudes. 
The  data  mode  for  each  force  run  Is  Identified  In  the  Test  Run  Summary  (Table 
2). 


Model  shadowgraph  or  schlieren  photographs  were  obtained  on  selected 
configurations  at  selected  model  attitudes  and  test  conditions  during  the 
pressure  phase  and  on  all  configurations  at  all  test  conditions  and  selected 
model  attitudes  during  the  force  phase. 


The  force  and  moment  measurements  were  reduced  to  coefficient  form  using 
the  digitally  filtered  data  points  and  correcting  for  first  and  second  order 
balance  interaction  effects.  Vehicle  coefficients  were  also  corrected  for 
model  tare  weight  and  balance-sting  deflections.  Model  attitude  and  tunnel 
stilling  chamber  pressure  were  also  calculated  from  digitally  filtered  values, 


Vehicle  aerodynamic  force  and  moment  coefficients  are  presented  in  the 
body-  and  stability-axis  systems.  Pitching  and  yawing  moment  coefficients  are 
referenced  to  a  point  on  the  model  centerline  24.0  inches  from  the  nose.  The 
stability-axis  system  coefficients  (CLS  and  CDS)  were  calculated  using  the 
forebody  axial  force  coefficient  (CA).  Model  diameter  and  base  area  were  used 
as  the  reference  length  and  area  for  the  aerodynamic  coefficients.  Model 
reference  dimensions  are  given  in  the  Nomenclature. 

3. 2. 2. 2  PWT  4T 

All  steady-state  measurements  were  sequentially  recorded 
by  the  facility  on-line  computer  system  and  reduced  to  the  desired  tinal  form. 
The  data  were  then  tabulated  in  the  Tunnel  4T  control  room,  recorded  on 
magnetic  tape,  and  transmitted  to  the  AEDC  central  computer  file.  The  data 
stored  in  the  central  computer  file  were  generally  available  for  plotting  and 
analysis  on  the  PWT  Interactive  Graphics  System  within  30  seconds  after  data 
acquisition.  The  immediate  availability  of  the  tabulated  data  permitted 
continual  on-line  monitoring  of  the  test  results. 

Surface  and  base  pressure  data  were  normalized  by  the  free  stream  static 
pressure  and  the  surface  and  pipe  pressure  data  were  reduced  to  pressure 
coefficient  form.  Selected  surface  pressure  data  were  also  presented  graph¬ 
ically  by  constructing  three-dimensional  color  contour  plots  over  the  model 
shape. 

The  model  force  and  moment  data  were  reduced  to  coefficient  form  in  the 
body-  and  stability-axes  systems.  The  model  reference  area  is  given  in  the 
Nomenclature  and  the  reference  lengths  are  given  in  Table  1.  The  moment 
reference  point  is  shown  in  Figure  2.  The  stability-axis  system  coefficients 
(CLS  and  CDS)  were  calculated  using  the  forebody  axial  force  coefficient  (CA) 
and  the  normal  force  coefficient  (CN).  The  base  pressure  and  its  area  (given 
in  Nomenclature)  were  used  to  calculate  the  base  axial  force. 


3.3  Uncertainty  of  Measurements 


3.3.1  VKF 

In  general.  Instrumentation  calibration  and  data  uncertainty 
estimates  were  made  using  methods  recognized  by  the  National  Bureau  of  Stan¬ 
dards  (NBS)  presented  In  Reference  10.  Measurement  uncertainty  Is  a  com¬ 
bination  of  bias  and  precision  errors  defined  as: 

U  *  ±(B  +  tg5  S) 

where  B  Is  the  bias  limit,  S  is  the  sample  standard  deviation,  and  tg5  Is  the 
95th  percentile  point  for  the  two  tailed  student's  "t"  distribution  (95 
percent  confidence  Interval),  which  for  sample  sizes  greater  than  30  is  taken 
equal  to  2. 

With  the  exception  of  the  force  and  moment  balance,  data  uncertainties 
are  determined  from  in-place  calibrations  through  the  data  recording  system 
and  data  reduction  program.  Static  load  hangings  on  the  balance  simulate  the 
range  of  loads  and  center-of-pressure  locations  anticipated  during  the  test, 
and  measurement  errors  are  based  on  differences  between  applied  loads  and 
corresponding  values  calculated  from  the  balance  equations  used  in  the  data 
reduction.  Lead  hangings  to  verify  the  balance  calibration  are  made  in  place 
on  the  assembled  model. 

Propagation  of  the  bias  and  precision  errors  of  measured  data  through  the 
calculated  data  were  made  In  accordance  with  Reference  10.  Uncertainties  for 
the  calculated  data  are  calculated  for  the  largest  measured  value  at  the 
primary  test  condition  on  each  parameter  at  each  Mach  number. 

3.3.2  4T 

The  aircraft  angles-of-attack  and  sideslip  were  corrected  for 
sting  deflections  caused  by  aerodynamic  loads.  The  flow  angularity  (AFA)  in 
the  tunnel  pitch  plane  was  determined  by  testing  the  aircraft  model  upright 
and  Inverted  and  the  flow  angularity  corrections  were  then  applied  to  the 
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data.  Corrections  for  the  components  of  model  weight,  normally  termed  static 
tares,  were  also  accounted  for  before  the  measured  loads  were  reduced  to 
coefficient  form. 


Uncertainties  (combinations  of  system  and  random  errors)  of  the  basic 
tunnel  parameters,  shown  in  Figure  6,  were  estimated  from  repeat  calibrations 
of  the  Instrumentation  and  from  the  repeatability  and  uniformity  of  the  test 
section  flow  during  tunnel  calibration.  Uncertainties  In  the  instrumentation 
systems  were  estimated  from  repeat  calibration  of  the  systems  against  secon¬ 
dary  standards  whose  uncertainties  are  traceable  to  the  National  Bureau  of 
Standards  calibration  equipment.  The  tunnel  parameter  and  Instrument 
uncertainties,  for  a  95  percent  confidence  level,  were  combined  using  the 
Taylor  series  method  of  error  propagation  described  in  Reference  10  to  deter¬ 
mine  the  uncertainties  of  the  parameters  In  Table  4. 


SECTION  4.0 


METHODOLOGY  DESCRIPTION 

4.0  Supersonic/Hypersonic  Arbitrary  Body  Program  (S/HABP) 

This  digital  computer  program  was  written  by  the  Douglas  Aircraft  Company 
and  is  documented  in  AFFDL-TR-73-159  (Reference  11).  The  program  is  a  com¬ 
bination  of  techniques  and  capabilities  necessary  in  performing  a  complete 
aerodynamic  analysis  of  supersonic  and  hypersonic  shapes.  The  program  was 
originally  designed  primarily  for  computing  the  aerodynamic  characteristics  of 
high-speed  arbitrary  reentry  vehicles.  Because  the  program  will  provide 
aerodynamic  coefficients  for  any  complex  arbitrary  shape  at  any  angle  of 
attack  for  Mach  numbers  equal  to  one  or  greater,  the  code  is  finding  renewed 
interest.  Many  of  the  new  state-of-the-art  missile  designs  are  now  non- 
circular/non-conventional  shaped  configurations  operating  at  higher  angles 
of  attack  and  supersonic  Mach  numbers. 

The  program  calculates  the  inviscid  aerodynamic  coefficients  by  using 
simple  pressure  coefficient  methods,  such  as  Newtonian  Impact,  on  the  external 
geometry  of  a  configuration  which  is  described  to  the  program  as  panels  of 
flat  quadrilateral  elements  as  shown  in  Figure  7a.  The  majority  of  the 
pressure  methods  are  simply  a  function  of  the  local  angle  of  attack  or  slope 
of  the  element  and  the  freestream  Mach  number.  The  viscous  or  skin  friction 
effects  are  computed  on  a  simplified  flat  plate  geometry.  Figure  7b,  using 
methods  such  as  Reference  Temperature.  The  program  has  15  different 
compression  (impact  angle  greater  than  0)  and  9  expansion  (impact  angle  less 
than  0)  methods  for  computing  inviscid  pressure  coefficients  plus  9  choices  of 
combinations  of  skin  friction  methods.  The  individual  forces  and  moments  of 
each  element  are  summed  up  by  the  program  to  give  component,  such  as  wing,  or 
tail,  aerodynamic  characteristics.  These  results  may  in  turn  be  summed  up  to 
give  parametric  buildup  or  complete  configuration  aerodynamic  characteristics. 
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4.2  Pressure  Integration 

An  integration  of  the  pressure  measurements  over  the  surface  of  the  model 
was  performed  to  determine  the  inviscid  aerodynamic  forces  and  moments.  The 
integrated  forces  and  moments  can  be  used  to  compare  directly  with  the 
inviscid  theoretical  calculations  and  by  subtracting  them  from  the  test  total 
force  and  moment  results  provide  the  increments  due  to  viscous  effects.  The 
instrumented  stations  did  not  extend  forward  of  station  3.2  so  a  theoretical 
value  of  the  nose  aerodynamics  was  generated  using  Tangent  Cone  and  Prandtl- 
Meyer  methods.  The  nose  analysis  used  the  S/HABP  code  with  the  methods 
described  above.  Between  station  3.2  and  36.0  the  pressure  data  was  integrated 
using  a  scheme  which  used  the  pressures  around  each  cross  section  to  determine 
a  value  at  each  station,  then  the  integration  proceeded  from  the  first  station 
to  the  base  to  produce  a  value  of  the  surface  integral  (Reference  12). 

The  pressure  at  each  tap  has  a  component  in  each  of  the  three  cartesian 
directions.  The  component  is  determined  from  the  direction  cosines  of  the 
surface  at  the  tap  location.  The  direction  cosines  were  calculated  from  the 
model  geometry.  The  direction  cosines  are  the  coefficients  of  the  unit  vector 
normal  to  the  surface.  The  unit  normal  vector  was  determined  by  the  vector 
product  of  two  orthogonal  vectors  defined  by  the  surface  geometry  functions. 

The  maximum  span  at  any  station  is  defined  as: 

a  =  amax  (x)1^ 

6 

and  the  maximum  height  at  any  station  is: 

b  =  bmax  (x)^ 

6 


The  cross  section  at  any  station  is  an  ellipse,  with  a  and  b  the  semima/jor  and 
semiminor  dimension.  The  remainder  of  the  dimensions  and  angles  may  now  be 
determined.  Referring  to  Figure  7c,  the  direction  cosines  of  the  surface  at  a 


point  P  are  evaluated  as  follows,  using  the  values  of  a  and  b  from  the  re¬ 
lations  shown  above  and  a  specified  value  of  the  meridian  angle  U. 

y  =  ((a2  *  b2)/(b2  +  a2tan2U))1/2 

z  =  y  *  tan  U 

W  =  atan  ( (b2/z)  *  (y/a2)) 

V  =  atan  (0.5  *  (y2  +  z2)1/2/x) 
the  two  orthogonal  vectors  are  therefore 

A  =  A^i  =  A2j  +  Agk 


B  =  B-j  i  +  B2j  +  Bjk 


A]  =  0 


A2  =  cos  W 


A3  =  -sin  W 


B1  =  cos  V 


B2  =  sin  V  *  cos  (90-U) 
B3  =  sin  V  *  sin  (90-U) 


The  unit  normal  at  the  surface  is  written  as 


"C  =  CjT  +  c2T  +  cjjfr 


and  the  direction  cosines  are 


C1  =  B2*  A3'  B3*A2 


The  force  and  moment  coefficients  may  be  determined  by  integrating  the  pres¬ 
sure  coefficients  over  the  surface.  The  formulation  for  the  longitudinal 
coefficients  is: 

CA  =  1/Ar  ff  Cp  *  Cj  ds  dh 

CN  =  1/Ar  ff  Cp  *  C3  ds  dh 

Cm  =  1/(Ar  *  1R)//  ((Cp*C1*2)-(Cp  *  C3  *  x))  ds  dh 

The  computerized  integration  routine  used  the  trapezoidal  rule,  first  around 
the  cross  section  on  the  distance  s,  then  down  the  body  along  the  length  h. 

Approximations  were  used  for  s  and  h.  The  chord  length  between  pressure 
taps  was  computed  as  the  arc  length.  They  were  summed  at  each  station  from 
the  top  center  to  the  bottom  center.  The  distance  h  was  computed  as  the  slant 
height  of  a  right  conic  solid,  with  s(x^)  and  s(x2)  being  the  periphery  of  the 
ends  of  the  solid. 

hl-2  =  ^x2  "  xl^  +  (s (x2 )/2ir  -  s(x1)/2tt)2 

4.3  NSWC  Euler  Code 

In  addition  to  comparing  pressure,  force,  and  moment  predictions  de¬ 
termined  with  the  Supersonic/Hypersonic  Arbitrary  Body  Program  to  the  experi¬ 
mental  data,  we  made  similar  comparisons  with  predictions  from  a  numerical 
computation  technique.  Only  inviscid  computer  codes  were  considered  since,  at 
the  time,  it  seemed  unwarranted  for  the  conditions  of  interest  (i.e.,  Mach 
number;  Reynolds  number;  and  body  configuration)  to  introduce  the  additional 
complications  of  a  viscous  method  such  as,  for  example,  a  Parabolized  Navier 
Stokes  Code  for  what  would  be  small  effects. 
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Experience  had  been  gained  on  one  particular  inviscid  code  during  previ¬ 
ous  in-house  studies.  This  code,  herein  called  the  NSWC  code,  is  a  forward 
marching  solution  to  the  steady,  inviscid,  supersonic  flow  equations.  It  was 
originally  written  at  the  Naval  Surface  Weapons  Center  (hence,  the  NSWC 
identification)  and  which  they  called  the  D3CSS  computer  code  (References  13 
and  14). 

We  need  an  initial  flow  field  data  plane  where  the  flow  field  is 
everywhere  supersonic  in  order  to  run  the  NSWC  code  for  a  body  in  supersonic 
flow.  For  the  original  use,  this  was  supplied  through  three  pre-programs 
called  BNT,  DDD,  and  BETA  which,  essentially,  provided  a  blunt  body  solution 
for  a  spherical  nose  at  angle-of-attack  (Reference  15).  The  complete  package 
of  codes  was  run  successfully  for  data  comparison  with  experimental  biconic 
data  and  the  results  reported  by  Scaggs  (Reference  16).  Use  of  the  blunt  body 
pre-programs  placed  a  lower  limit  on  free  stream  Mach  number  of  about  4  and, 
since  the  Mach  number  range  of  the  data  on  the  elliptical  power  law  bodies  was 
between  1.76  and  5.03  an  alternative  solution  was  needed. 

The  code  chosen  to  provide  the  starting  solution  was  the  CM3DT  computer 
program,  written  by  Science  Applications  International  Corporation  (SAIC)  for 
the  Ballistic  Missile  Office  (BMO)  (Reference  17).  The  code  is  a  time- 
dependent,  steady-state  solution  for  supersonic/hypersonic  flow  over  nosetips 
of  arbitrary  shape  and  yields  the  asymptotic  limit  of  the  unsteady  flow 
problem.  The  NSWC  code  was  modified  to  accept  input  from  the  CM3DT  code  by 
SAIC  as  part  of  an  earlier  contract  effort  (Reference  18). 

For  all  cases  produced  for  this  report,  a  perfect  gas  condition  was  used 
and  calculations  were  stopped  at  an  angle  of  attack  of  10  degrees  since  the 
NSWC  code,  like  other  supersonic  inviscid  codes,  fails  if  either  the  axial 
velocity  component  becomes  subsonic  or  if  axial  flow  separation  occurs.  This 
seemed  to  be  at  just  above  10-degrees  angle  of  attack,  especially  for  the 
higher  ellipticity  body,  since  the  12-degree  angle-of-attack  case  always 
failed.  Finally,  although  the  bodies  of  interest  were  symmetrical  ellipses  in 
cross  section,  this  particular  option  did  not  exist  in  the  code.  Therefore, 
the  bi-ellipse  option  was  used  with  the  0=0°  and  180°  values  of  y  the  same 


» 


at  each  X  location.  Ten  planes  of  body  geometry  were  required  to  properly 
describe  the  nose  shape  up  to  an  axial  location  1  inch  from  the  nosetip. 

4.4  FL057  Euler  Code 

FL057  is  a  finite  volume  Euler  method  which  was  modified  to  permit 
arbitrary  geometries  through  the  use  of  multiple  grid  blocks.  The  volume  or 
cells  are  defined  by  eight  neighboring  grid  points.  Conservation  of  mass, 
momentum,  and  energy  are  satisfied  in  an  integral  form  on  each  volume  by  a 
pseudo-time-stepping  four-stage  Runge-Kutta  scheme.  Flow  variables  are 
assumed  to  be  located  at  cell  centers,  permitting  centered  differences  to  be 
used  to  provide  second-order-accurate  spatial  derivatives.  Artificial  dis¬ 
sipation  is  added  to  suppress  the  odd-even  point  decoupling  that  is  typical  of 
center-differenced  Euler  methods  and  also  to  reduce  non-physical  pressure 
oscillation  around  shocks  and  stagnation  points.  The  dissipative  terms  are 
calculated  by  blending  fourth  and  second  differences  and  are  scaled  by  second 
derivatives  of  pressure.  To  increase  convergence  rates,  both  enthalpy  damping 
and  implicit  residual  smoothing  is  used.  Surface  boundary  condition  is 
normal-flow  imposed  using  only  the  cell  adjacent  to  the  surface.  Far  field 
boundary  conditions  are  of  a  nonreflecting  type  (Reference  19).  The  FL057  is 
a  time-dependent  solution  that  allows  solutions  to  be  obtained  at  subsonic  as 
well  as  supersonic  speeds. 

We  did  not  attempt  to  generate  the  missile  grid  to  align  the  3-D  grid 
with  the  bow  shock  shape  since  the  grid  was  used  for  a  range  of  Mach  numbers 
and  angles  of  attack.  The  bow  shock  is  dependent  on  the  configuration  angle 
of  attack  and  freestream  Mach  number.  Shock  smearing  will  occur  when  the  bow 
shock  is  unaligned  with  the  grid,  introducing  an  unknown  amount  of  error  into 
the  solution.  At  the  nose  of  the  configuration  the  shock  would  approximate 
the  shape  of  the  blunt  nose  and  therefore  align  with  the  grid  at  moderate 
supersonic  Mach  numbers  (Reference  20). 

4.5  Missile  Datcom 

Missile  Datcom  provides  an  aerodynamic  design  tool  with  the  predictive 
accuracy  suitable  for  preliminary  design,  yet  has  the  utility  to  be  extended 


by  rapid  substitution  of  methods  to  fit  specific  applications.  The  code  uses 
a  component  build-up  approach  to  calculate  the  static  stability  and  control 
characteristics  of  missiles  with  both  unconventional  fin  arrangements  and 
arbitrary  cross  sections.  The  primary  advantage  of  component  build-up  methods 
over  panel  methods  Is  speed  of  operation. 
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SECTION  5.0 


TEST  RESULTS 

5.1  VKF  Tunnel  A 

5.1.1  Force  and  Moment  Testing 

The  effect  of  ellipticlty  ratio  on  the  basic  aerodynamic  charac¬ 
teristics  is  shown  in  Figures  8  and  9  for  the  elliptical  bodies  at  Mach  2.0 
and  5.0.  The  trends  in  the  normal  force  and  pitching  moment  coefficients  are 
consistent.  Increasing  the  ellipticity  of  the  body  increases  the  normal  force 
coefficient,  lift-to-drag  ratio  and  the  pitching  moment  coefficient.  At 
angles  of  attack  less  than  5  degrees  the  axial  force  coefficient  increases 
with  increasing  ellipticity  ratio.  Interestingly,  at  higher  angles  of  attack, 
the  2.5:1  body  shows  a  higher  axial  force  coefficient  than  either  the  2.0:1  or 
the  3.0:1.  In  fact,  in  most  cases,  the  3.0:1  ratio  body  axial  force  coeffi¬ 
cient  approaches  the  value  of  the  2.0:1  body  at  higher  angles  of  attack. 

Figures  10  and  11  show  the  effects  of  Mach  number  on  the  3.0:1  body  at  a 
length  Reynolds  number  of  two  million  per  foot.  For  clarity,  five  Mach 
numbers  are  plotted  in  each  set  of  figures;  Mach  2.0,  2.5,  3.0,  4.0,  and  5.0 
in  one  set  and  Mach  1.76,  2.0,  2.5,  3.0,  and  3.5  in  the  other.  The  plots  show 
that  the  normal  force  and  pitching  moment  coefficients  decrease  with  increas¬ 
ing  Mach  number.  At  small  angles  of  attack  (less  than  6  degrees)  axial  force 
coefficient  decreases  for  increasing  Mach  number.  The  axial  force  coefficient 
increases,  however,  for  increasing  Mach  number  at  higher  angles  of  attack. 

This  corresponds  to  a  decrease  in  lift-to-drag  ratio.  This  trend  reversal 
occurs  at  higher  angles  of  attack  for  the  2.0:1  and  2.5:1  bodies,  and  is 
attributed  to  the  formation  of  large  leeside  vortices. 

The  effects  of  Reynolds  number  are  shown  in  Figures  12  and  13.  As 
expected,  this  range  of  change  in  Reynolds  number  had  virtually  no  effect  on 
either  normal  force  or  pitching  moment  coefficients  at  either  Mach  2.0  or  5.0. 
The  effects  of  Reynolds  number  on  axial  force  coefficient  varied  with  Mach 
number.  At  Mach  2.0  (Figure  12),  the  axial  force  at  zero  angle  of  attack  was 


only  slightly  higher  at  the  lower  Reynolds  number,  but  the  difference  in¬ 
creases  with  increasing  angle  of  attack.  This  is  surprising  since  the  only 
difference  in  axial  force  should  be  due  to  skin  friction  which  should  be 
relatively  constant.  At  Mach  5.0  (Figure  13),  axial  force  coefficient  de¬ 
creases  as  expected  for  increasing  Reynolds  number  at  low  angles  of  attack. 
Above  10-degrees  angle  of  attack,  axial  force  coefficient  for  the  3.0:1  body 
is  essentially  independent  of  Reynolds  number. 

5.1.2  Pressure  Results 

During  the  theoretical/experimental  comparison  effort,  plots  were 
made  of  pressure  ratio  vs  x-station  for  the  top  and  bottom  centerlines  at  zero 
angles  of  attack  for  all  the  test  Mach  numbers.  These  plots  showed  some 
unexpected  trends  in  the  pressure  data.  At  zero  angle  of  attack,  one  would 
expect  symmetrical  pressure  distribution  for  the  top  and  bottom  centerlines. 
The  data,  however,  showed  jumps  in  the  pressure  coefficient  midway  down  the 
body  on  both  the  top  and  bottom  (Figure  14).  Consultation  with  the  test 
engineers  at  AEDC  confirmed  suspicions  that  the  fluctuations  were  tunnel 
induced,  and  could  be  expected  for  any  pressure  test  at  very  low  angles  of 
attack.  Further  investigation  showed  that  the  variations  in  small  values  of 
pressure  have  a  negligible  effect  when  measuring  and  integrating  pressures  at 
larger  angles  of  attack. 

Figure  15  shows  representative  plots  of  pressure  coefficient  versus  body 
radial  angle  for  all  11  x-stations.  Data  shown  is  for  3.0:1  body  at  12- 
degrees  angle  of  attack  at  Mach  2.0  and  5.0. 

5.1.3  Flow  Visual ization 

5. 1.3.1  Vapor  Screens 

Results  of  the  vapor  screen  test  are  shown  in  Figure  16 
for  the  3:1  ellipticity  model  at  Mach  3  and  16-degrees  angle  of  attack.  Both 
upstream  and  downstream  views  are  included.  Figure  17  includes  upstream  views 
for  12-  and  20-degrees  angle  of  attack  as  well  (Reference  6).  These  composite 
photos  show  the  location  and  growth  of  the  shocks  and  leeside  vortices.  The 
darker  areas  in  the  photographs  indicate  less  water  vapor 


than  the  lighter  regions.  Several  explanations  for  this  phenomenon  have  been 
proposed.  One  explanation  is  that  the  high  rotational  velocities  of  the 
vortex  ejects  the  water  droplets.  Another  is  that  the  water  vapor  does  not 
penetrate  through  the  shear  layer  Into  the  boundary  layer,  which  feeds  the 
vortex.  Another  possibility  is  that  the  water  undergoes  a  phase  change  and 
vaporizes  due  to  the  temperature  increase  produced  by  the  shock  (Reference 
21).  A  good  description  of  the  basic  techniques  used  for  vapor  screen  testing 
is  contained  in  Reference  22. 

For  the  large  angles  of  attack  shown  in  the  figures,  the  shocks  are 
clearly  visible  with  the  leeside  vortices.  Secondary  vortices  are  also 
visible,  especially  in  the  upstream  views.  They  are  located  just  inboard  of 
the  leading  edge,  under  the  sheet  feeding  the  primary  vortex. 

5. 1.3. 2  Shadowgraph/Schlieren 

Samples  of  the  shadowgraph  and  schlieren  data  are  shown 
in  Figures  18  and  19.  Both  are  for  the  3:1  ellipticity  model  at  12-degrees 
angle  of  attack.  The  shadowgraph  data  is  Mach  3.0  and  the  schlieren  photo  is 
for  a  Mach  2.5  run.  Figure  20  shows  the  comparisons  of  the  experimentally 
obtained  shock  shape  with  the  NSWC  Euler  code.  The  code  predicts  the  com¬ 
pression  shock  very  well;  on  the  expansion  side,  however,  the  code  underpre¬ 
dicts  the  shock  angle  especially  at  higher  angles  of  attack. 

5. 1.3.3  Oil  Flows 

Samples  of  the  oil  flow  photographs  are  shown  in  Figure 
21  for  the  3:1  configuration  at  Mach  3.0  and  12-degrees  angle  of  attack. 

Figure  22  shows  the  primary  and  secondary  separation  angles,  measured  graph¬ 
ically  from  the  photographs.  The  primary  separation  angle  is  defined  as  the 
line  of  separation  just  past  the  leading  edge  as  the  flow  moves  from  the 
windward  to  the  leeside  of  the  body.  The  secondary  separation  angle  is  the 
separation  line  further  inboard  on  the  leeside  where  the  flow  reattaches. 


Figure  22  shows  that  for  all  angles  of  attack,  the  primary  separation 
point  gets  closer  to  the  leading  edge  as  the  flow  moves  down  the  body.  In 
general,  the  radial  separation  angle  increases  with  increasing  angle  of 
attack.  The  trends  are  reversed  for  secondary  separation  angles.  The  angle 
gradually  decreases  as  the  flow  moves  down  the  model,  and  separation  angle 
decreases  for  increasing  angle  of  attack.  The  data  drops  sharply  in  the  nose 
region,  probably  due  to  the  initial  formation  of  the  vortices.  The  complex 
flow  in  the  nose  region  is  evident  in  the  oil  flow  photographs. 

5.2  PWT  4T 

5.2.1  Force  and  Moment  Test 

The  effects  of  ellipticity  ratio  on  the  longitudinal  aerodynamic 
characteristics  are  shown  in  Figures  23  through  25  for  Mach  numbers  0.4,  0.8, 
and  1.3. 

The  trends  in  the  test  data  are  consistent  for  varying  ellipticity  ratio 
and  Mach  number.  Increasing  the  ellipticity  of  the  body  increases  normal 
force  coefficient,  1 ift-to-drag,  and  the  pitching  moment  coefficient.  The 
axial  force  coefficient  generally  increases  slightly  with  increasing  elli¬ 
pticity  ratio,  the  exception  being  at  higher  angles  of  attack  and  at  the 
lowest  Mach  numbers,  where  it  shows  either  a  very  small  decrease  or  no  change 
at  all.  The  lateral  directional  derivatives  are  also  consistent  for  changing 
ellipticity  ratio  (Figure  26).  The  rolling  moment  derivative  becomes  increas¬ 
ingly  negative  with  increasing  ellipticity  ratio  and  the  yawing  moment  deriva¬ 
tive  becomes  less  negative  at  all  Mach  numbers.  The  side  force  coefficient  is 
more  positive  with  increasing  ellipticity  ratio  at  all  angles  of  attack  for 
Mach  0.4  and  0.55.  This  trend  Is  reversed  above  10-degrees  angle  of  attack  at 
the  higher  Mach  numbers,  where  the  side  force  coefficient  decreases  with 
increasing  ellipticity  ratio. 

5.2.2  Pressure  Test 

The  transonic  tunnel  force  and  moment  data  showed  dramatic  changes 
in  side  force  and  yawing  moment  with  varying  angle  of  attack  below  Mach  0.6 
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and  4-degrees  yaw  angle.  The  side  force  coefficient  goes  from  -0.02  to  +0.06 
between  8-  and  14-degrees  angle  of  attack,  while  the  yawing  moment  coefficient 
changes  slope  and  even  goes  positive  for  the  3.0:1  model  (Figure  26).  Analy¬ 
sis  of  the  pressure  coefficient  data  showed  large  changes  in  pressure  dis¬ 
tribution  with  changing  angle  of  attack,  particularly  in  the  negative  pressure 
coefficients  about  the  leading  edge.  To  illustrate  this.  Figure  27  shows  a 
polar  plot  of  Cp  versus  body  radial  angle  for  the  3.0:1  body  at  6-  and  12- 
degrees  angle  of  attack  and  Mach  0.4  at  the  16-inch  x-station.  The  plot  shows 
that  at  6  degrees  the  negative  pressure  coefficients  balance  out  and  a  nega¬ 
tive  side  force  results  from  the  unbalanced  positive  pressure  coefficients  on 
the  bottom  surface.  At  12  degrees  the  positive  pressure  coefficients  are 
slightly  unbalanced,  but  the  negative  Cp  on  the  one  side  is  much  higher  than 
the  other,  resulting  in  a  positive  side  force.  These  same  types  of  changes 
caused  by  the  vortices  on  the  upper  surface  result  in  a  change  in  the 
lengthwise  pressure  distribution  which  affects  the  yawing  moment  as  well. 


SECTION  6.0 


DATA/ PREDICT I ON  COMPARISONS 
6.0  Force  and  Moment  Comparisons 

6.1.1  Supersonic/Hypersonic  Arbitrary  Body  Program  (S/HABP) 

Some  typical  elliptical  body  S/HABP  theoretical  results  compared 
with  the  test  data  are  shown  in  Figures  28  and  29,  (3.0:1  body  at  Mach  numbers 
2.0  and  5.0).  These  theoretical  results  use  the  S/HABP  Tangent  Cone  pressure 
method  for  compression  surfaces  and  Van  Dyke  Unified  method  for  expansion. 
Turbulent  skin  friction  was  selected  for  the  viscous  computations.  At  lower 
Mach  numbers  the  only  coefficient  predicted  with  reasonable  accuracy  is 
lift-to-drag  ratio.  The  overprediction  of  both  normal  force  and  axial  force 
compensated  each  other  to  give  a  reasonable  value  for  L/D.  Above  Mach  3.0  the 
agreement  between  normal  force,  pitching  moment,  and  L/D  was  fairly  good. 

Axial  force  coefficient  was  still  slightly  overpredicted,  however. 

6.1.2  Transonic  Force  and  Moment  Data 

Typical  Missile  Datcom  theoretical  results  compared  with  the 
static  longitudinal  force  and  moment  test  data  are  shown  in  Figures  30  through 
32  for  the  3.0:1  ellipticity  ratio  configuration  at  Mach  numbers  0.4,  0.8,  and 
1.3.  These  theoretical  results  used  the  second-order  shock-expansion  method 
for  axisymmetric  bodies  at  supersonic  speeds  and  turbulent  skin  friction  at 
all  Mach  numbers.  The  Missile  Datcom  did  a  good  job  of  predicting  normal 
force  coefficient  below  10-degrees  angle  of  attack,  where  CN  is  fairly 
linear.  At  higher  angles  of  attack,  it  underpredicted  the  values.  At 
subsonic  Mach  numbers  the  predicted  pitching  moment  coefficient  is  about  twice 
the  test  data  value  at  corresponding  values  of  normal  force  coefficient.  At 
Mach  1.3  the  predicted  pitching  moment  values  agree  very  well  with  test  data 
even  to  the  highest  angles  of  attack.  Because  the  base  axial  force 
coefficient  is  such  a  large  portion  of  the  body  total  axial  force  coefficient, 
both  of  these  coefficients  have  been  plotted  and  compared  with  predicted 
values.  The  predicted  axial  force  coefficients  due  to  the  base  pressure  are 


low  at  all  three  Mach  numbers,  particularly  at  Mach  1.3.  The  total  axial 
force  coefficient  is  fairly  well  predicted  at  Mach  0.4,  but  is  low  at  Mach  0.8 
and  very  low  at  Mach  1.3  where  it  drops  off  drastically  with  increasing  angle 
of  attack.  The  lift-to-drag  ratio  predictions  are  not  too  bad  at  Mach  numbers 
of  0.4  and  0.8,  but  are  high  at  Mach  1.3. 

Comparisons  of  the  Missile  Datcom  with  the  lateral  directional  coeffi¬ 
cients  for  the  3.0:1  ellipticity  ratio  configurations  at  Mach  numbers  0.4, 

0.8,  and  1.3  are  shown  in  Figures  33-35.  The  present  version  of  Missile 
Datcom  gave  no  values  for  rolling  moment  coefficient.  The  side  force 
coefficient  was  predicted  very  well  at  low  angles  of  attack  for  all  three  Mach 
numbers.  The  yawing  moment  was  predicted  as  being  more  negative  in  value  than 
the  test  data  at  Mach  0.4  and  0.8  and  slightly  less  negative  at  Mach  1.3.  The 
test  data  shows  a  reversal  in  sign  for  the  side  force  coefficient  near  8 
degrees  angle  of  attack,  going  from  negative  values  of  -  0.01  to  positive 
values  as  high  as  0.06. 

6.1.3  FL057  Euler  Code 

Figure  36a  presents  a  comparison  between  measured  and  predicted 
values  of  normal  force  and  pitching  moment  coefficient  versus  angle  of  attack 
for  the  2.5:1  configuration  at  Mach  0.55.  Below  6-degrees  angle  of  attack 
FL057  does  an  excellent  job  of  predicting  CN  and  CM.  Above  6-degrees,  where 
the  data  is  no  longer  linear,  the  code  underpredicts  the  test  data.  The 
differences  between  test  and  prediction  in  both  plots  is  evidence  that  the 
influence  of  vortices  is  not  present  in  the  Euler  results. 

Figure  36b  shows  comparisons  for  the  same  configuration  at  Mach  2.0.  The 
predicted  and  test  values  are  in  excellent  agreement  below  6- degrees  angle  of 
attack.  Above  6-degrees,  the  slopes  are  the  same  but  are  shifted  by 
approximately  0.5-degrees  angle  of  attack.  One  possible  explanation  for  this 
is  the  smearing  of  the  bow  shock  and  grid  (Reference  20).  The  versus  Cm 
curve  throws  excellent  agreement  throughout  the  angle-of-attack  range. 
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6.2  Integrated  Pressures 

Figures  37  and  38  show  the  comparison  of  the  integrated  pressure  forces 
and  moments  with  the  force  and  moment  test  data.  Also  plotted  are  the  theo¬ 
retical  results  from  the  NSWC  Euler  code.  The  agreement  between  the  integrat¬ 
ed  normal  force  and  pitching  moment  coefficients  and  the  force  and  moment  test 
data  was  excellent.  The  forces  and  moments  on  the  nose  to  the  X  *  3.2  station 
were  calculated  with  the  S/HABP  code  using  the  Tangent  Cone/Prandtl -Meyer 
methods,  and  were  then  added  to  the  test  pressure  forces  and  moments.  The 
excellent  agreement  indicates  that  the  theoretical  methods  used  for  the  nose 
section  are  quite  good.  The  Increment  in  the  axial  force  coefficient  between 
the  integrated  pressure  data  and  the  force  test  data  is  due  to  the  viscous 
forces;  i.e.,  the  skin  friction.  The  plots  show  almost  a  constant  increment 
in  axial  force  coefficient  due  to  skin  friction  until  about  10-degrees  angle 
of  attack  where  the  increment  starts  to  increase.  The  increase  in  the  axial 
force  coefficient  increment  is  very  large  at  the  highest  angles  of  attack  for 
the  lower  Mach  numbers.  As  an  example,  the  increment  of  axial  force 
coefficient  goes  from  about  0.05  at  8-degrees  angle  of  attack  to  0.11  at  20- 
degrees  angle  of  attack  for  the  3-to-l  ellipticity  ratio  configuration  at  Mach 
2.0.  Figure  39  shows  the  axial  force  coefficient  increment  for  the  3.0:1 
ellipticity  ratio  configuration  for  Mach  numbers  2.0,  3.0,  4.0,  and  5.0.  Also 
shown  on  the  figures  are  the  skin  friction  predictions  from  the  S/HABP  code. 
The  S/HABP  skin  friction  methods  calculate  about  the  right  value  for  the  skin 
friction  increment  at  low  angles  of  attack,  but  the  increment  is  essentially 
constant  with  increasing  angle  of  attack.  At  the  highest  angles  of  attack  for 
all  three  configurations  the  predicted  axial  force  coefficients  due  to  skin 
friction  were  less  than  half  what  the  test  data  indicated. 
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6.3  Cp  versus  Body  Radial  Angle 


6.3.1  Supersonic/Hypersonic  Arbitrary  Body  Program  (S/HABP) 


Typical  S/HABP  theoretical  pressure  results  compared  with  test 
data  at  x-stations  of  3.2,  16.0,  and  35.2  inches  are  shown  in  Figures  40  and 
41  at  Mach  2.0  and  5.0  for  the  3:1  body.  The  S/HABP  theoretical  results  shown 
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are  for  the  Tangent  Cone  method  for  compression  surfaces  and  Van  Dyke  Unified 
method  for  expansion. 

The  plots  show  that  the  bottom  surface  pressures  are  generally  predicted 
very  well  by  the  Tangent  Cone  method  until  the  leading  edge  of  the  body  is 
approached;  there  the  test  pressure  coefficients  drop  off  and  become  negative 
around  10  degrees  before  the  leading  edge  at  0  =  90  degrees.  The  S/HABP  code 
calculates  positive  local  deflection  angles  for  all  points  on  the  lower 
surface  of  the  body  up  to  and  including  the  leading  edge  when  the  body  is  at 
positive  angle  of  attack  (Figure  42).  None  of  the  compression  methods  in  the 
S/HABP  code  which  are  functions  of  local  deflection  angle  and  Mach  number  will 
give  negative  pressure  coefficients  for  positive  deflection  angles.  The  plots 
show  that  the  test  pressure  coefficients  at  angle  of  attack  have  fairly  high 
negative  values  on  the  top  in  the  vicinity  of  the  leading  edge  which  are  more 
negative  than  predicted.  The  predicted  values  further  up  on  the  top  reach  the 
maximum  negative  values  of  the  test  and  remain  constant  across  the  rest  of  the 
upper  surface  to  the  center  of  the  body,  while  the  test  data  becomes  much  less 
negative,  even  slightly  positive  in  some  cases,  on  the  leeside  away  from  the 
leading  edge. 

The  large  negative  pressure  coefficients  about  the  leading  edges  of  the 
elliptical  bodies  strongly  influences  their  aerodynamic  characteristics.  The 
importance  of  the  leading  edge  pressures  is  directly  related  to  the  amount  of 
surface  area  over  which  they  act.  The  ratio  of  wetted  area  ( A^)  to  total  band 
area  (Ay)  for  a  0.2-inch-length  band  of  body  about  X  =  16  inches  is  plotted 
versus  the  body  angle  in  Figure  43  for  the  2-to-l  and  3-to-l  ellipticity 
ratio  bodies.  The  plot  shows  that  for  a  10-degree  increment  in  body  angle  the 
amount  of  wetted  area  near  the  leading  is  proportionally  larger  so  that  the 
large  negative  pressures  near  the  leading  edge  (o  =  80  to  100  degrees)  are 
acting  over  a  very  larye  surface  area. 


TJWJV.^ 

r. 


•J 


6.3.2  NSWC  Euler  Code 

Comparisons  of  the  pressure  coefficients  predicted  by  the  NSWC 
Euler  code  for  x-stations  of  3.2,  16.0,  and  35.2  are  also  shown  in  Figures  40 
and  41  for  which  the  program  provided  values.  The  plots  show  that  the  NSWC 
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program  does  a  very  good  job  of  predicting  the  pressure  coefficients  about  the 
body.  The  program  predicts  both  the  negative  pressure  coefficients  on  the 
compression  side  just  before  the  leading  edge  and  the  drop  off  in  negative 
pressure  coefficients  on  the  leeside  away  from  the  leading  edge  wall.  The 
good  agreement  with  test  data  indicates  the  potential  for  these  types  of  codes 
in  predicting  elliptical  body  aerodynamics.  Techniques  for  extending  their 


range  of  applicability  and  methods  for  calculating  viscous  forces  for  inclu¬ 
sion  with  the  invicid  results  should  be  investigated  further. 

6.3.3  FL057  Euler  Code 


Figure  44  shows  comparisons  of  pressure  coefficient  versus  span 
for  three  x  stations;  3.2,  16.0,  and  35.2.  Data  shown  are  for  the  2.5:1  model 
at  Mach  2.0  and  12-degrees  angle  of  attack.  At  the  x=3.2  station  the  shapes 
of  the  curves  are  the  same  but  the  FL057  prediction  is  shifted  in  the  negative 
Cp  direction.  At  the  x=16.0  station,  the  code  still  predicts  an  attached  flow 
condition,  while  the  flow  visualization  test  data  shows  the  formation  of 
vortices  on  the  upper  surface.  This  effect  is  also  present  at  the  x=35.2 
station.  The  shift  in  Cp  may  be  due  to  a  smeared  bow  shock,  which  creates 
angularity  in  the  flow  and  results  in  an  apparent  angle-of-attack  change 
(Reference  20). 


6.4  Cp  versus  Local  Deflection  Angle 

Since  the  prediction  methods  in  the  S/HABP  code  are  functions  of  the 
freestream  Mach  number  and  the  local  deflection  angle,  a  set  of  charts  was 


prepared  to  show  the  relation  between  the  measured  pressure  coefficient  and 
the  calculated  local  deflection  angle.  The  charts  show  the  extent  of  corre¬ 
lation  with  local  deflection  angle  and  reveal  the  difficulty  of  devising  a  new 
pressure  function  which  would  predict  the  aerodynamics  more  accurately. 


The  pressure  coefficients  are  shown  with  the  local  deflection  angle  in 
Figures  45  and  46.  Each  graph  shows  all  the  measured  values  on  the  model  at  a 
single  angle  of  attack.  Representative  data  was  selected  at  angles-of-attack 
increments  of  approximately  4  degrees.  The  data  symbol  indicates  the  general 
region  of  the  pressure  tap  as  being  on  the  top,  the  leading  edge,  or  the 


bottom  of  the  model.  Tangent  Cone  theory  is  also  marked  on  the  graphs.  The 
highest  values  of  the  pressure  coefficient,  which  occur  on  the  bottom  center, 
are  very  near  the  theory.  At  locations  other  than  the  bottom  center  the 
pressure  coefficient  is  less  than  theory.  For  these  tests,  the  Tangent  Cone 
theory  represents  an  upper  bound. 

The  effect  of  Mach  number  is  most  evident  at  negative  surface  deflection 
angles.  The  measurements  on  the  top  of  the  model  indicate  a  flow  field  which 
is  influenced  by  parameters  other  than  the  local  angle.  The  range  of  pres¬ 
sures  on  the  top  indicate  the  effects  of  boundary  layer  separations  and  the 
development  of  vortices  in  the  flow  field.  At  Mach  2.0  the  pressures  are  very 

sensitive  to  boundary  layer  separation.  The  pressure  coefficients  have  values 

2 

from  Cp  =  0  to  a  lower  bound  near  Cp  =  -1/M  .  The  lower  bound  is  based  on 
past  work  using  base  pressure  measurements  from  wind  tunnel  and  flight  tests 
which  has  shown  that  the  maximum  attainable  suction  pressure  is  about  7-tenths 

vacuum.  The  equation  for  pressure  coefficient  as  pressure  goes  to  zero  is 

2  2 
-2/a  M  which,  for  a  =  1.4,  when  multiplied  by  0.7  equals  -1/M  . 

The  correlation  with  local  deflection  angle  also  shows  that  the  leading 
edge  pressures  are  much  lower  than  Tangent  Cone  theory  would  predict.  The 
difference  is  larger  at  low  Mach  numbers,  but  still  significant  at  Mach  5. 


SECTION  7.0 

RESULTS  AND  CONCLUSIONS 

An  extensive  data  base  for  elliptical  cross  section  bodies  has  been 
generated  for  use  in  missile  design  activities  as  well  as  establishing  a 
benchmark  for  the  evaluation  of  aerodynamic  performance  prediction  programs. 
The  comprehensive  combination  of  force,  pressure  and  flow  visualization  data 
will  allow  identification  of  the  source  of  deficiencies  in  current  analysis 
techniques  and  indicate  improvements  to  be  made.  This  data  base  has  been 
documented  in  a  series  of  in-house  Technical  Memorandums. 

During  this  study  a  number  of  different  types  of  analysis  codes  have  been 
used  to  generate  theoretical  results  for  data/theory  comparisons.  Both  of  the 
Euler  codes  used  in  this  study,  FL057  and  NSWC,  did  a  very  good  job  of  pre¬ 
dicting  the  pressure  distribution,  normal  force  and  pitching  moment  at  angles 
of  attack  of  6  degrees  or  less.  Since  the  Euler  codes  are  inviscid  the  axial 
forces  predicted  do  not  include  the  viscous  effects  and  the  correct  prediction 
of  the  axial  force  coefficient  requires  that  the  viscous  effects  be  accounted 
for.  The  large  differences  between  the  integrated  test  pressure  axial  force 
coefficients  and  total  force  axial  force  coefficients  show  the  importance  of 
the  viscous  forces,  particularly  at  higher  angles  of  attack.  Analysis  methods 
for  predicting  viscous  effects  other  than  just  simple  strip  theory  skin 
friction  calculations,  as  in  S/HABP,  will  be  required.  Evaluations  of  the 
importance  of  the  effects  of  flow  field  vortices,  boundary  layer  separation, 
and  boundary  layer  transition  should  be  done.  Computational  Fluid  Dynamics 
(CFD)  codes  such  as  a  Parabolized  Navier-Stokes  code  need  to  be  evaluated  for 
this  class  of  configuration. 
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Figure  1.  Tbnnel  A 
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TABLE  1.  Model  Configuration  Designation 


CONFIGURATION 


DESCRIPTION 


B20 


2.0:1  ELLIPTICAL  BODY,  amax  =  4.162  in. 


bmax  —  2.081  in. 
L  =  8.324  in. 


B25 

2.5:1  ELLIPTICAL  BODY,  amax  =  4.654  in. 

bmax  =  1-862  in. 

R! 

L  =  9.308  in. 

gh 

B30 


3.0:1  ELLIPTICAL  BODY,  amax  =  5.098  in. 


br^ax  =  1699  in. 
L  =  10.196  in. 


37 


$ 
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TABLE  2.  TEST  RUN  SUMMARY 
VKF  Hinnel  A  Run  Log 


FITCH  •t’N  AT  CONSTANT  BE^A 


TABLE  2.  TEST  RUN  SUMMARY  (Continued) 

VKF  Hinnel  A  Run  Log 
Fbrce  and  Moment  Test 


M 

ALPHA 

*E  6 
*10“6 

PITCH  tfS  AT  CONSTANT 

_ °  i  2  i  4 

3.51 

A1 

2.0 

68 

A2 

69 

4.02 

A1 

2.0 

44 

A2 

45 

46 

4.51 

A1 

2.0 

47 

A2 

48 

5.03 

A1 

2.0 

63 

A2 

! 

64 

65 

5.04 

A1 

3.0 

66 

A1 

4.0 

67 

NOTES: 


1.  tins  for  BETA  •  0  wr*  run  la  continuous  sweep  aode,  except  as  noted, 
and  tor  BETA  •  2  and  4  vara  run  la  point-pause  aode. 

2.  *  Indicates  point-pause  aode. 


3.  ALPHA  schedule:  A1  -  -4, -3. 5, -3,-2. 5, -2, -1.5, -1,-0. 5, 0,0. 5, 1,1. 5, 2, 2. 5, 

3,3.5,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20 
A2  -  -4,-2,-1,0,1,2,4,6,8,10.12,16.20 
A3  -  -2, -1.5. -1,-0. 5, 0,0. 5, 1,1. 5, 2, 2. 5. 3, 3. 5, 4, 5, 6, 7, 8, 
9,10,11,12,13,14,15,16,17,18,19 


.y. 


TABLE  2.  TEST  RUN  SUMMARY  (Continued) 
VKF  Tbnnel  A  Run  Log 

Shadovgraph/Schlleren  Flovflcld  Photographic  Log 


PDL  Elliptic  Bodies  Test 


BB/rr  boll 


CODE 

BUN 

M 

«10* 

•0. 

ALPHA 

P.S. 

30 

2 

3.01 

3.0 

496 

-0.18 

4.32 

8.66 

13.01 

17.53 

21  88 

1-6 

3 

2.0 

-0.23 

4.19 

8.39 

12.61 

16.95 

21.15 

6 

1.0 

-0.20 

4.08 

8.17 

12.26 

16.46 

20.54 

7 

2.50 

2.0 

505 

•0.21 

4.23 

8.50 

12.80 

17.21 

21.47 

12 

2.00 

1.0 

496 

-0.22 

4.12 

8.26 

12.44 

16.76 

20.91 

13 

1.76 

2.0 

465 

-0.13 

4.28 

8.60 

13.01 

17.61 

1-5 

25 

14 

1.76 

2.0 

-0.14 

4.22 

8.48 

12.81 

17.34 

1-5 

15 

2.00 

1.0 

-0.18 

4.10 

8.20 

12.35 

16.6? 

20.74 

1-6 

16 

2.0 

-0.23 

4.20 

8.47 

12.79 

17.28 

21.61 

19 

2.50 

2.0 

-0.20 

4.18 

8.42 

12.68 

17.07 

21.31 

21 

3.01 

3.0 

718 

-0.25 

4.26 

8.55 

12.87 

17.34 

21.64 

22 

2.0 

-0.24 

4.15 

8.34 

12.52 

168  6 

21.02 

25 

1.0 

-0.24 

4.05 

8.13 

12.20 

16.41 

20.46 

20 

26 

3.01 

1.0 

-0.22 

4.05 

8.11 

12.18 

16.36 

20.40 

27 

2.0 

-0.21 

4.10 

8.25 

12.40 

16.69 

20.81 

30 

3.0 

-0.19 

4.19 

8.42 

12.66 

17.06 

21.29 

31 

2.50 

2.0 

455 

-0.23 

4.16 

8.34 

12.53 

16.88 

21.07 

36 

2.00 

1.0 

718 

-0.19 

4.07 

8.16 

12.25 

16.50 

20.60 

37 

1.76 

2.0 

-0.13 

4.17 

8.38 

12.62 

17.09 

1-5 

3B 

4.02 

2.0 

-0.27 

4.04 

8.13 

12.21 

16.44 

20.52 

1-6 

25 

41 

4.02 

2.0 

-0.18 

4:09 

8.19 

12.31 

16.55 

20.65 

30 

44 

4.02 

2.0 

-0.19 

4.11 

8.22 

12.35 

16.62 

20.73 

47 

4.51 

2.0 

-0.18 

4.08 

8.18 

12.28 

16.52 

20.61 

25 

49 

4.51 

2.0 

-0.18 

4.07 

8.15 

12.23 

16.45 

20.53 

■>r» 

•»  •» 

51 

4.51 

2.0 

-0.19 

4.06 

8.14 

12.20 

16.41 

20.45 

53 

5.03 

2.0 

-0.16 

4.03 

8.09 

12.16 

16.34 

20.38 

56 

5.04 

3.0 

-0.16 

4.07 

8.17 

12.27 

16.51 

20.61 

57 

4.0 

-0.18 

4.12 

8.25 

12.41 

16.71 

20.86 

25 

59 

5.04 

3.0 

-0.17 

4.10 

8.22 

12.35 

16.64 

20.76 

r  •*» 

W  w 

5.03 

2.0 

8.15 

16.40 

20.47 

1-4 

30 

63 

5.03 

2.0 

-0.18 

4.07 

8.15 

12.24 

16.47 

20.54 

1-6 

66 

5.04 

3.0 

-0.16 

4.12 

8.25 

12.39 

16.68 

20.83 

67 

4.0 

-0.14 

4.19 

8.39 

12.58 

16.95 

21.17 

68 

3.51 

2.0 

-0.23 

4.14 

8.33 

12.51 

16.85 

21.02 

25 

70 

3.51 

2.0 

-0.21 

4.13 

8.27 

12.43 

16.74 

20.88 

20 

72 

3.51 

2.0 

-0.20 

4.09 

8.22 

12.35 

16.62 

20.74 

MOTES:  1.  BUN  Mo.  and  Photo  Sequence  Mo.  (P.S.)  are  first  and  second  number 
respectively,  on  photographs. 
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TABLE  2.  TEST  RUN  SUMMARY  (Continued) 
VKF  Ttmnel  A  Pressure  Tfest 


ALPHA 


wa  *wVv>o£h 


N  N  »n  «  H  <« 
OA  H  N  n  lA  <0 


H  >6  H 
00  O'  pH 

00  00  on 


<#  tA  «0  CD 


SJA  O  O  n  O  sA  IA  ON  N  to  00  N  <  00  H  N  N  no  ON  <N  »A  on  <  00  H  »A  O'  N  iA  00 

O'  H  N  n  tA  »o  H  ON  HN  <  in  O  CO  N  00  O  fN  Os  ©  pH  0  «0  m  N  n  lA  sA  fN  O  pH 

«•  O'  ON  O'  O'  ON  CM  pH  pH  pH  pH  pH  pH  pH  ©  O  pH  O'  O'  O  O  O  O  O 


n  n  ia 
P4  N  N 


NO  N  O'  O  H 
N  N  N  H  CA 


oo  n  co  n  no  O'  cs  oO  -o 

N  ONO  N  NO  o  O' 

00  00  O'  O'  O'  O'  O'  CM  pH 


®  H  is 
HN-O 
HH  HHj 


lA  ®  H  ^  N  CM  fN  © 

SpH  cn  »n  n  n  m 

O  O  O  O  CM  CM  CM 


00H-ON 
NO  N  On  Q  H 
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TABLE  2.  TEST  RUN  SUMMARY  (Continued)  , 

VKF  Tunnel  A  -  Run  log 

Vapor-Scrten  Test 
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TABLE  2.  TEST  RUN  SUMMARY  (Continued) 

VKF  Tunnel  A  -  Run  log 

Phase  1  Oil-Flow  Test 


TABLE  2.  TEST  RUN  SUMMARY  (Continued) 

VKF  Tunnel  A  -  Run  log 

Phase  2  Oil  Plow  Test 


MACH  ALFnA  OIL  VISCOSITY 

RUN  NO.  MODEL  deg  TOP/BOTTOM 
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3 

it 
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10 

11 

12 


3.0  B30 

< 1 
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i 1 
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12 

14 
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14 

12 

10 
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Light. 'Mediun 
Heavy /Heavy 
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14 
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10 
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18  2.0 

19 

20 
21 
22 
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10 
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23  I 

24  t 


B30 
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14 

12 
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Light  oil  was  80X  10  Centl  Stokes  (CS)  and  201  100CS 
Medium  oil  was  501  IOCS  and  501  100CS 
Heavy  oil  was  1001  1000CS 


ALPHA  Schedule:  AU  -4.-3,-?.-! .0,1 3, 4, 6,8, 10, 1?, 16,20  deg 
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CONFIG  M  BETA 


b.  Radial  Location 

Pressure  Orifice  Locations  (Continued) 
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Table  J.  Pressure  Orifice  Location  and  Designation 


3.2 
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12.8 

16.0 

19.2 
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25.6 

28.8 

32.0 

35.2 
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14 

27 

40 

59 

78 

97 

116 

135 

154 

173 

20 
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15 

28 

41 

60 

79 

98 

117 

136 
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30 

42 

61 

80 

99 

118 

137 

156 
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40 
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16 

29 

43 

62 

81 

100 

119 

138 

157 

176 

SO 

44 

63 

82 

101 

120 

139 

158 

177 

55 

a 

17 

30 

45 

64 

83 

102 

121 

140 
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60 

46 

65 

84 

103 

122 
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179 

65 
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18 

31 

47 

66 

85 
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123 

142 
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70 

48 

67 

86 
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143 
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75 

6 

19 

32 

49 

68 

87 
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125 

144 
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182 

00 

69 

88 

107 

126 
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85 

7 

20 

33 

51 

70 

89 
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127 

146 
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184 

90 

52 

71 

90 

109 

128 

147 
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185 

95 

8 

21 

34 

53 

72 

91 

110 

129 

148 
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105 

9 

22 

35 

54 

73 

92 
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149 

168 

187 

115 

10 

23 

36 

55 

74 

93 

112 

131 

150 

169 

■ 

188 

125 

11 

24 

37 

56 

113 

132 

151 
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189 

140 

12 

25 

38 

57 

76 

95 

114 

133 

152 

171 

190 

ISO 

13 

l6 

39 

58 

77 

_ 

96 

_ 

115 

134 

153 

— 

172 

191 

NOTE:  Thermocouples  located  at  approximately  X  «  22.4  and  THETA  «  0  and  180. 
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Figure  5.  TUnnel  Model  Installation 


b.  Tlinnel  4T 

Figure  5.  Tlinnel  Model  Install 
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ELLIPTICITY  RATIO  2*1 


a.  Inviscid  Geometry 
Figure  7.  S/HABP  Geometry 


Pressure  Integration  Direction  Cosines 
ure  7.  S/HABP  Geometry  (Concluded) 
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Figure  II.  Mach  Number  Effects  (Concluded) 
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Figure  12.  Reynolds  Number  Effects,  M  =  2.0  (Concluded) 
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Figure  15.  Pressure  Coefficient  About  Body 
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Figure  15.  Pressure  Coefficient  About  Body  (Concluded) 
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Figure  16b.  VAPOR  SCREEN  PHOTO 


Figure  19  Schlieren  photograph 
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Figure  24.  Ellipticity  Ratio  Effects,  M  =  0.8  (Concluded) 
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Figure  25.  Ellipticity  Ratio  Effects,  M  =  1.3  (Concluded) 
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Figure  26.  Stability  Derivatives,  M  =  0.4 


ELLIPTICAL  BODY  MISSILE 

3.0:1  ELLIPTICITY  RRTIO 
MflCH  2.00 

RUN  9  TEST  DRTR 
S/HRBP  THEORY  -  5,5 


RTTRCK  RNGLE  OF  RTTRCK 


ELLIPTICAL  BODY  MISSILE 

3.0:1  ELLIPT1CITY  R FIT  10 
MflCH  5.03 


RTTRCK 


ELLIPTICAL  BODY  MISSILE 

3.0:1  ELLIPTICITY  RATIO 
ARCH  5.03 


ELLIPTICAL  BODY  MISSILE 

MACH  0.40 

REYNOLDS  NUMBER  -  2.40X10+6  PER  FT 
3.0=1  ELLIPTICITY  RATIO 


-10.0  0.0  10. 0 

ANGLE  OF  ATTACK  ANGLE  OF  ATTACK 


ELLIPTICAL  BODY  MISSILE 

MACH  0.80 

REYNOLDS  NUMBER  -  2.39X10+6  PER  FT 
3.0:1  ELLIPTIC  I TY  RATIO 


h->  CQ 

cr  cc 
o  o 

o  < 


<3  O 
O 


ro  g  < 

Q<3 


OO'I-  CT- 

°  O 


figure  31.  Fbrce  And  Moment  Comparisons,  M  =  0.8  (Concluded) 
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Figure  32.  Fbrce  And  Moment j Corripa rison s,  M  =  1.3  (Concluded) 
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Figure  37.  Integrated  Pressure  Comparisons,  M  =  2.0  (Concluded) 
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Figure  38.  Integrated  Pressure  Comparisons,  M  =  5.0  (Continued) 
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Figure  38.  Integrated  Pressure  Comparisons,  M  =  5.0  (Concluded) 
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Figure  40.  Cp  vs.  Phi  Angle  Comparisons,  M  =  2.0  (Continued) 
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Figure  40.  Cp  vs.  Phi  Angle  Comparisons,  M  =  2.0  (Continued) 
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Figure  40.  Cp  vs.  Phi  Angle  Comparisons,  M  =  2.0  (Continued) 
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Figure  40.  Q>  vs.  Phi  Angle  Comparisons,  M  =  2.0  (Concluded) 
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no 


VvSV 


v.v 


rngle: 


Figure  41.  Cp  vs.  Phi  Angle  Comparisons.  M 


Z'D 

IJ30 


133 


PRESSURr  rill  it.  PRESSURE  COOT  I'KO.iSlJRr;  COOT, 

o.o  0.2  o.i  -o  , »  0.0  0.2  n.*  -0.2  O.'.l  0.2 


SURFACE  LOCAL  DEFLECTION  ANGLE 


Data 


1.17 


PRESSURE  COEFFICIENT 


CP  VS.  LOCRL  DEFLECTION  RNGLE 


3.0:1  ELLIPTICITY  RRTIO 
MRCH  NO.  -  2.00 


RUN  NO.  -  143 
MODEL  rlphr  -  0.0 


D  -  BOTTOM 
O  -  LERDING  EDGE 
A 


TANGENT  CONE  THEORY 


P*  Ptri 


prhssurf:  cofjticient 


PRESSURE  COEFFICIENT 


/vv.v; 


CP  VS.  LOCHL  deflection  rngle 

3.0:1  ELLIPTICITY  RRTIO  RUN  NO.  -  148 
MRCh  NO.  -  2.0G  MODEL  At-PHR  -  8.4 


□  -  BOTTOM 
C  -  LEADING  EDGE 
A  "  TOP 


TANGENT  CONE  THEORY 


PRESSURE 


CP  VS.  LOCAL  deflection  angl 

3.0:1  ELLIPTICITY  RRTIC 

RJN  NO.  -  1=0 

MRCH  NO.  -  2.00 

CD 

O 

MODEL  RLFHR  -  12.8 

D  -  BOTTOM  d- 

i 

TANGENT  CONE  THEORY 

3  -  LERDING  EDGE 

1 


CP  VS.  LOCAL  DEFLECTION  ANGLE 

3.0:1  ELLIPTICITY  RATIO  RUN  NO.  -  152 
MACH  NO.  -  2.00  MODEL  AlPHA  -  17.2 

CD 

o'-| 


D  -  BOTTOM 
0  -  LEADING  EDGE 


u) 

b 


TANGENT  CONE  THEORY 


CP  VS.  LOCRL  DEFLECTION  RNGL 

3.0:1  ElL!p7  4  '—  I  T  -  RATIO  RUN  NO.  -  154 

MPO.  NC.  -  2. GO  MODEL  Rl=^A  -  21.6 

© 

o  " 


3  -  BOTTOM 
3  -  LEADING  EDGE 


ID 

o 


TANGENT  CONE  THEORY 


PRESSURE  COEFFICIENT 


CP  vs.  local  deflection  angle 


3.0:1  ELLIPTICITY  RRTIO 
MRCH  NO.  -  5.03 


RUN  NO.  -  255 
MODEL  RlPHR  -  0.0 


□  -  BOTTOM 
O  -  LERDING  EDGE 
A  -  TOP 


TANGENT  CONE  THEORY 


-24.0 

-1/M2 


CP  VS.  LOCRL  DEFLECTION  RNGLG 

3.0:1  ELLIPTICITY  RRTIO  RUN  NO.  -  25S 
MRCH  NO.  -  5.03  MODEL  RLPNR  -  4.1 


LOCRL  RNGLE  OP  RTTRCK 


Figure  46.  Cp  vs.  Local  Deflection  Angle,  M  =  5.0  (Continued) 
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Figure  46.  Cp  vs.  Local  Deflection  Angle,  M  =  5.0  (Concluded) 
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NOMENCLATURE 


A,  Ar 
AE 

ALPHA 

d 

a 

max 

BETA 

b 

b 

max 

CA,  CA 

CAB 

CAT 

CDS,  CQ 

CLL ,  C1 

CLM,  Cm 

CLN ,  Cp 
CLS ,  CL 
CN,  C, 
CODE 
CONFIG 
Cp 

CY,  Cy 
L,  1R 
L/Ds 
LM 
M 

MU 


2 

Reference  area,  27.210  inch 

Base  area,  27.210  inch2 

Angle  of  attack,  degree 

Semimajor  (horizontal)  span  at  X,  inch 

Semimajor  span  at  model  base,  in.  (see  Table  1) 

Sideslip  angle,  degree 

Semiminor  (vertical)  height  at  X,  inch 

Semiminor  height  at  model  base,  inch  (see  Table  1) 

Forebody  axial  force  coefficient,  body  axes,  CAT-CAB 

Base  axial  force  coefficient,  body  axes,  -(PBA-P)AB/Q*A 

Total  axial  force  coefficient,  body  axes,  total  axial  force/Q*A 

Drag  coefficient  (based  on  CA),  stability  axes 

Rolling  moment  coefficient,  body  axes,  rolling  moment/Q*A*L 

Pitching  moment  coefficient,  body  axes,  pitching  moment,  Q*A*L 

Yawing  moment  coefficient,  body  axes,  yawing  moment/Q*A*L 

Lift  Coefficient  (based  on  CA),  stability  axes 

Normal  force  coefficient,  body  axes,  normal  force/Q*A 

Model  configuration  number 

Model  configuration  designation 

Pressure  coefficient,  (P-P  /Q) 

oo 

Side  force  coefficient,  body  axes,  side  force/Q*A 
Reference  length.  Diameter  (inch) 

Lift-to-drag  ratio  (based  on  CA),  stability  axes 
Model  length,  36.000  inch 
Free-stream  Mach  number 

9 

Dynamic  viscosity  based  on  free-stream  temperature,  lbf-sec/ft‘~ 


NOMENCLATURE  (Continued) 


NCP 

P 

PBA 

PB  i 

psfa 

psia 

psid 

PHI 

PREF 

PT 

PT2 

Q 

RE 

RHO 

RUN 

T 

TT 

V 

X 

XMRP 

YCP 


Normal  force  center-of-pressure  location,  body  axes,  inches  from 
nose;  XMRP-(CLM-L/CN)  or  XMRP-(CLM-AO*L/CN-AO)  for  ALPHA  =  0  and 
BETA  =  0 

Free-stream  static  pressure,  psia 

Average  base  pressure  (PBT  +  PBB  +  PBL  +  PRR)/4,  psia 

Base  pressure,  i  =  T,  B,  L,  and  R,  where  T,  B,  L,  and  R  are  top, 
bottom,  left,  and  right  looking  upstream,  respectively,  psia 
Pound  per  square  foot  area 

Pound  per  square  inch  area 

Pound  per  square  inch  differential 

Roll  angle,  degree 

Reference  pressure  for  ZOC  pressure  module,  psia 

Tunnel  stilling  chamber  pressure,  psia 

Total  pressure  downstream  of  a  normal  shock,  psia 

Free-stream  dynamic  pressure,  psia 

Free-stream  unit  Reynolds  number,  ft"* 

3 

Free-stream  density,  Ibm/ft 

Data  set  identification  number 

Free-stream  static  temperature,  °R 

Tunnel  stilling  chamber  temperature,  °R,  or  °F 

Free-stream  velocity,  ft/sec 

Axial  location  from  nose  of  model,  inch 

Axial  distance  from  model  nose  to  model  moment-reference  lo¬ 
cation,  24.000  inch 

Side  force  center-of-pressure  location,  body  axes,  inches  from 
nose,  XMRP  -  (CLN*L/CY) 
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